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We study the late-time behavior of a universe in the framework of brane gas 
cosmology. We investigate the evolution of a universe with a gas of super- 
gravity particles and a gas of branes. Considering the case when different 
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1 Introduction 



Unifying gravity with other forces of nature strongly suggests that there may 
be more than three spatial dimensions in the unifying scale. Some of these 
are hidden from the low energy observers. This is related to the cosmological 
question why we have only three spatial dimensions. There had been earlier 
attempts to answer this question |TJ El El HI E] • One of the notable is the idea 
of string cosmology proposed by Brandenberger and Vafa (BV) This is a 
mechanism to generate dynamically the spatial dimensionality of spacetime 
and to explain the problem of initial singularity. The key ingredient of this 
model is based on the symmetry of string theory called T-duality. With 
this symmetry the spacetime has a topology of nine dimensional torus and 
its dynamics is driven by a gas of fundamental strings. Cosmology based 
on this string-modified Einstein-Friedman equations was studied in many 
directions El El EH • 

Developments in string theory during the last decade revealed that it is 
not a theory of only strings and has richer structure of branes jTU] . Recently 
there have been attempts to understand the cosmological evolution in this 
new framework of string theory with D-branes E| and the mechanism 
of BV was considered in [HI]. This model of brane gas cosmology was stud- 
ied extensively [HI [T31 HH1 E] . In the picture of brane gas cosmology, the 
universe starts from hot, dense gas of D-branes in thermal equilibrium. The 
winding modes of branes obstruct the growing of the spatial dimensions. 
Branes with opposite winding numbers can annihilate if their world volume 
interacts. Thus hierarchy of scales can be achieved between the wrapped 
and unwrapped dimensions. Of particular interest is whether the unwrapped 
configuration of branes can successfully inflate to make the unwrapped di- 
mensions grow indefinitely while making the wrapped ones remain small or 
at least grow much slowly [111 Ej • 

In the authors studied the late-time cosmology in M-theory with a 
supergravity particle gas and wrapped 2-brane gas. In this paper, we study 
the late-time behavior of brane gas cosmology in the context of string theory 
by extending the formalism of ^B] to general p-brane gas. We investigate the 
behavior of a universe with a gas of supergravity particle and a gas of wrapped 
p-brane. Considering the case where different directions are anisotropically 
wrapped by various branes, we will argue the possible hierarchical evolution 
of scales between wrapped and unwrapped dimensions. 
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2 Brane gas dynamics 

In this section, we extend the formalism of brane gas dynamics in jT^j for 
branes of arbitrary dimensionality and set up some preliminaries for our cal- 
culation. The model we will consider is type II string theory compactified on 
T 9 . Let us consider the late stage of BV scenario where the radii and curva- 
ture scales are grown larger than the ten-dimensional Planck length. When 
the radii are grown enough, we can neglect the brane-antibrane annihilation. 
And the excitations on the branes will be red-shifted away faster than the 
brane tension. Supergravity is a good approximation with the growing radii 
and falling temperature. Matter fields can be classified into two types. One 
is the massless supergravity particles corresponding to bosonic and fermionic 
degrees of freedom. We ignore massive modes since these will decay quickly. 
The other is wrapped D-branes. These branes do not interact with each other 
when the radii are grown large. We also ignore the fluctuations within the 
brane. We assume that the supergravity particles and brane gases are homo- 
geneous for simplicity of calculation. Then we can take the averages of the 
contributions of all kinds of particles and branes to the energy-momentum 
tensor. 

In the point of string theory, the gravitational interaction is described by 
the coupled system of the metric and dilaton. Dilaton plays an important role 
in the large-small symmetry of string theory called T-duality. To be specific, 
let us start from the following effective action of type II string theory 



where <p is the dilaton field and S m denote all matter actions including the 
brane action. Since we are interested in the late time cosmology, we work in 
the Einstein frame, defined by 



We drop the superscript E from now on. The equations of motion from this 
action are 




(1) 



9mn — e2 9mn- 
In terms of the Einstein metric, the action can be written as 



(2) 




(3) 
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(5) 



f=g 5(f) 

Since the string coupling (g = e^) is considered to be small in our assumption 
of large radii, we do not consider the running of the dilaton throughout this 
paper. Then the equations of motion are described simply by the Einstein 
equation 



R 



1 



MN 



9mnR 



1 5S„ 



2 — - ^g- 6g MN- 

With the metric ansatz of torus, with D — 9, 

D 

-dt 2 + ^2(oi(t)) 2 del 0<9i< 2tt, 
i=i 



(6) 



(7) 



the non-vanishing components of the Einstein tensor are 




where i is not summed in the second equation. 

For the matter part, we first consider a gas of massless supergravity par- 
ticles, with energy density ps an d pressure p$. Since we assume the gas to be 
homogeneous and isotropic, we take a perfect fluid form of energy momentum 
tensor 



T jv = diag(-ps,Ps,---,Ps)- (10) 

The equation of state appropriate for D spatial dimensions fixes ps = (1/D)ps- 
The second source of energy momentum comes from a gas of branes, 
wrapped on the various cycles of the torus. The matter contribution of a 
single p-brane to the action, in string frame, is represented by the Dirac- 
Born-Infeld (DBI) action 
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Sp — T n 



V p J d p+1 ^^-det( giIU + + 2ira>F» v ), (11) 
where T p is the tension of p-brane and is the induced metric to the brane 

dX M Q X N 

g^u - Qmn d ^ ■ (12) 

Here M, N are the indices of (D + 1) dimensional bulk spacetime and fi, v 
are those of brane. B^ u is the induced antisymmetric tensor field and 
is the field strength tensor of gauge fields living on the brane. With the 
assumption of grown radii, we neglect B^ u and F^ u terms below. Ignoring 
the running of the dilaton, we can absorb the effect of constant dilaton into 
the redefinition of brane tension in the Einstein frame. 

Let us consider a gas of 1-branes. The contribution of these branes are 
characterized by wrapping numbers iVj and iVj, where we take iVj to represent 
the number of 1-branes wrapped on the % cycle, while Ni represents the 
number of anti- 1-branes. By symmetry it is a reasonable assumption that 
they are equal. The wrapping numbers are based on the understanding 
of thermal fluctuations in the early universe. A single 1-brane action is 
described by 

Sl-brane = ~ T l j ^a/ - det 9 a/3, (13) 

where T\ is the 1-brane tension. The stress tensor from this 1-brane action 
is 

t mn = _ Ti J d 2£ 5 D+i( X _ X(0) V 7 " det g g af3 d a X M dpX N . (14) 

As an example, for a single 1-brane wrapped on the 1 cycle and uniformly 
smeared over the transverse T " 1 , the stress energy tensor is 



T M N = — diag(l, 1, 0, 0, . . . , 0). (15) 

2710.2 ' ' ' 2najy 

With wrapping numbers Ni and iVj for i cycle, the non-zero components of 
the 1-brane gas stress tensor are 
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T\ = -y27rJ2 a k(N k + N k ), (16) 

k 

T\ = -^2na t (N t + N t ). (17) 

Now we insert these two sources of energy-momentum into the right hand 
side of the Einstein equations, 

G M N = —8nGT M N . (18) 

The time component equation can be solved for the energy density of the 
super gravity gas, 

8nG Ps = £ £ — " ^(2vr) £ a k {N k + N k ). (19) 
z T7, a k a i v , 

kjtl k 

The space component equations, after some algebra, can be reduced to the 
following set of second-order differential equations 



di 8-kGTx 



_D + 1 
D 



a t V 

1 d k di ^ dfcdj 

kytl kjtl 

For a gas of 2-branes, the wrapping can be characterized by JV^, where 
we take A^<j to represent the number of 2-branes wrapped on the (ij) cycle, 
while Niyj represents the number of anti- 2-branes. The stress energy tensor 
from a single 2-brane action is 



t mn = _ T2 J d S£ S D+i( X _ X(0) v 7 - det g g a(3 d a X M d p X N , (21) 

where T 2 is the 2-brane tension. Then, for a single 2-brane wrapped on the 
(12) cycle and uniformly smeared over the transverse T D ~ 2 , the stress energy 
tensor is 



With the wrapping numbers, the non-zero components of stress energy tensor 
from the 2-brane gas are 



T\ = -y(27i) 2 J2^aiN kl , (23) 
k^i 

T\ = -y(2n) 2 J2aka t (N kt + N lk ). (24) 
Repeating the same procedure as in 1-brane case, we have 

8nG Ps = J E — " E ( 25 ) 



hi 8nGT 2 
a< 1 V 



^ D + \ ^ E a k a i N ki - (2?r) 2 E afc °* ( Nki + 



k^l kj±i 

For a gas of 3-branes, we characterize the wrapping numbers of (ijk) cycle 
as Nij k , where we take with = 1 to represent the number of 3-branes 
wrapped on the (ijk) cycle, while Nij k with €ij k = —1 represents the number 
of anti-3-branes. The stress energy tensor for a single 3-brane action is 

t mn = _ Ts J d ^ D+ \X - X{0) y/- det g g af3 d a X M dpX N , (27) 

where T 3 is the 3-brane tension. For a single 3-brane wrapped on the (123) 
cycle and uniformly smeared over the transverse T D ~ 3 , the stress energy 
tensor is 

T M N = -z — diag(l, 1, 1, 1, 0, . . . , 0). (28) 

The non-zero components of the 3-brane gas stress energy tensor are 
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T t 



-^(2vr) 3 y - 

V K ' ^ 3 

k^l^m 



Q-kO'lO'mNklm, 



T\ = -y(2vr) 3 \ a kaia l {N ikl + N llk + 



Nkii) 



(29) 
(30) 



k^l^i 



where Ylk^i^m means that the sum is taken when all three indices are differ- 
ent. The constraint equation and the second-order differential equations for 
the radii are 

8nGps = ~ V] - 87T ^ Ts (2?r) 3 V] \a k aia m N k i m} (31) 
2 dfcOi V . 'rrf 3 



kj^l^m 



Oi 
CI: 



8ttCT 3 



3L> + 1 



— — — (2tt) 2^ -^akaia m N k i m 

k^l^m 



(27r) 3 ^ ^a^a* (JV iW + N hk + iV fcK ) 



2jD a ^ rr' a ^ a * 



(32) 



3 Late time behavior 

To study the solutions of the field equations (J2UJ) . (J2T?)) and (J3~2"j) . we introduce 
new variables Aj(t) = ln(27raj(t)) as in 0. The field equations when the 
matter fields are coming from supergravity particles and wrapped 1-brane 
gas become, from (J2UJ) . 



\ + = 8nG f^p s + ^jPb) - ^^e A <(iV 4 + #.), (33) 

where the volume of the spatial torus V can be expressed as V = Xi and 
Pb is the energy density from 1-brane tension 



Tl 



PB = y^e Afc (iV l + iV fc ). 



(34) 



The energy density from supergravity particles is fixed by the constraint 
equation (fTT?|) 

87rGps = ^A l A J -87rGp B . (35) 

This system of equations can be regarded as a non-relativistic particle 
moving in D dimensions. The particle has a coefficient of friction, given by 
V/V, due to the expansion of the universe. One can consider the right hand 
side of (JBlIj) gives two position-dependent forces acting on the particle [TE]. 
The first term from the supergravity particles, 

if } = SttG (Ipg + ^Pb) , (36) 

is positive definite. This force is common to every value of % and drives a 
uniform expansion of the universe. The second term from the brane gas 
contribution, 

F« = -^e^ + % (37) 

can be either zero or negative depending on the wrapping number. This term 
can suppress the growth of dimensions when wrapping numbers are nonzero. 
Thus anisotropic wrapping numbers can lead to an anisotropic expansion of 
the dimensions. 

The constraint equation ()35|) can be rewritten as 

I>) =^(A,) 2 + 167tG( Ps + pb). (38) 

i J i 

The right hand side is positive definite for all nontrivial cases and we choose 
the direction of time to make V > 0. 

We are interested in cases where the wrapping number Ni is anisotropic. 
We can classify the spatial dimensions into two kinds. We refer to a direction 
i as unwrapped if N = N = and as wrapped if Ni = N ^ 0. Now consider 
the case when m dimensions are unwrapped and D — m dimensions are 
wrapped. This configuration can be achieved by thermal fluctuations in the 
early universe through the mechanism of Brandenberger and Vafa [3] . Denote 
the two subvolumes as in JH] 
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7/1 



fj, = Xi = In (volume of unwrapped torus), (39) 
i=i 

D 

A = Xi = In (volume of wrapped torus). (40) 

j=m+l 

Summing over the appropriate values of i, and using the definition of p&, we 
find the following differential equations for p, and A 



: x „ / m £iu \ 

H+(ji + A)fi = 8nG f — p s + -jjZTi pB J ' ( 41 ) 

v / ana ^f D ~ m D-2m + l \ 

A + (A + A)A = SttG f — ps + ^ _ - PbJ- (42) 

Repeating the same procedure when the matter fields are coming from 
supergravity particles and wrapped 2-brane gas, we have, from (J26)) . 



\ + ^\ = 8nG ( ip s + ^Pb) - £(Jfy + Nji)e Xj . (43) 

The constraint equation has the same form as ()35|) with the energy density 
Pb replaced by 

PB = y$>\^AV (44) 

In this case we can classify the spatial dimensions into three classes jTHj. A 
direction i is referred to as unwrapped if Nj = Nji = for all j. Directions 
% for which and Nji are nonzero except for those j corresponding to an 
unwrapped direction are referred to as fully wrapped. Directions % where 
some of the Nj or Nji are zero for values of j which are not unwrapped are 
referred to as partially wrapped. The equations governing the motion of \i 
and A are 



ji + (fi + A)/i 
A + (/i + A) A 



.. r , . m 3m 



. ri ,'D-m D-3m + 2 
87rG I — - — p s H — Pb 



D 



D- 1 



(45) 
(46) 
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Similarly when the matter fields are coming from supergravity particles 
and wrapped 3-brane gas, we have, from (J22J), 



•■ V • / 1 4 

K + y K = 8irG I — p s + jj—[Pb 



8nGT 3 \ 
-e A 



\e^e Xk {N i]k + N kl] +N ]ki ). (47) 



The constraint equation has the same form as with the brane energy 
density ps replaced by 

PB = yYl \^e x *e x «N ijk . (48) 

Also we can classify the spatial dimensions into three classes as in 2-brane 
case; unwrapped, partially wrapped and fully wrapped. The equations for p 
and A are 



(m Am \ 

D ps + ~D^T pB J ' (49) 

A + (/i + A)A = 8ttG( — — p s + D _ 1 Pb) ■ (50) 

Generalizing the above result to the case when matter fields are coming 
from supergravity particles and wrapped p-brane gas, one can find 



p + (p + A)p = 8nG fe s + ^^Pb) ' (51) 

: s : ^ f D — m D — (p + l)m + p \ . . 
A + (p + A)A = SvrG^— ^Ps + ^{ -PbJ- (52) 

The evolution of unwrapped subvolume p and wrapped subvolume A depends 
on the number of unwrapped dimensionality m. For small m, both terms on 
the right hand side of A equation (|52~jl are positive. If A is zero, then the 
second derivative must be positive, leading to a local minimum. Conversely, 
for a larger value of m, the right hand side of this equation has both a 
positive and a negative term. In this case both local maxima and minima 
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are possible. However, for unwrapped subvolume jl, there are only positive 
terms on the right hand side. So if these directions are initially expanding 
they will expand forever. 

When the second term of right hand side in A equation vanishes the 
brane tension does not contribute to the growth of the internal dimensions. 
This gives a criterion for the critical spatial dimensionality of the unwrapped 
subspace 

D — 2m c +1 = 0, for 1 — brane, 

D-3m c + 2 = 0, for 2 - brane, (53) 

D — 4m c + 3 = 0, for 3 — brane. 

For D — 9, this occurs at m c = 5, 11/3 and 3 for 1-, 2- and 3-brane gases 
respectively. From this result we can conclude that wrapping the internal 
dimensions with 1- and 2-branes cannot yield the observed three large spatial 
dimensions. For 3-brane there is no contribution from p-Q but A will grow by 
ps- However the growing rate is smaller than that of p. When p is wrapped 
by p-dimensional brane gas, the criterion is 

D — (p + l)m c + p = 0, for p— brane. (54) 

As an example, for p — 4, this gives m c = 13/5. This means that when 
m > 13/5 there can be negative contribution from the p-Q terms. Our analysis 
implies that the compact extra dimensions should be wrapped with three or 
higher dimensional branes for the observed three large spatial dimensions. If 
they are wrapped by just one or two dimensional branes, they cannot remain 
compact. 

4 Conclusion and Discussion 

We have studied the late-time behavior of the universe in the framework of 
brane gas cosmology. For the cases when different dimensions are anisotrop- 
ically wrapped by various branes, we have derived Friedman-like equations 
governing the dynamics of wrapped and unwrapped subvolumes. We pointed 
out that one cannot keep the extra dimensions compact with one or two di- 
mensional brane gases. 

Though the wrapped directions grow slowly compared with the unwrapped 
dimensions, it is not guaranteed that the compact dimensions can be stabi- 
lized. If four or higher dimensional branes are frozen at the end of thermal 
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stage, these branes give negative force to A while giving positive force to 
fx. If we assume that every term in /?b takes a comparable contribution to 
the energy density, we can replace the effect of all wrapped branes with a 
(D — m)-brane gas. In this point of view, one can study the evolution of a 
universe starting from the metric 



where a 2 is the scale factor of the three dimensional space and b is the scale 
factor of the internal D — 3 dimensional subspace. The stabilization of the 
extra dimensions with this asymmetric setting was attempted in the brane 
world scenario ^H] ■ Recently the stabilization of extra dimensions was stud- 
ied numerically. In string gas cosmology, cosmological solution for the late 
time evolution showed that large dimensions continue to expand and the 
small ones are kept undetectably small \F3\. In (201 > ^ is shown that six 
compact dimensions become stabilized at the self-dual radius while three di- 
mensions grow large. More recently it has been shown that, by reducing the 
effect of string gas to the four dimensional Einstein gravity, string modes can- 
not stabilize the internal dimensions wrapped by winding strings, except in 
the special case of one extra dimension [23] . However, the evolution of extra 
dimensions is slow enough compared to the unwrapped large dimensions. 

One way to achieve the stability of the internal dimensions is to include 
the effect of fluxes. Alexander ^5] pointed out the possibility of stabilizing 
the extra dimensions with G-flux in the context of brane gases in eleven di- 
mensional supergravity derived from M-theory. We expect similar argument 
can be applied to brane gases in ten dimensional supergravity from type II 
string theory. For example, if we consider the effect of NS-NS sector in the 
bulk action, the cosmology is described by 



which involves the dynamics of two-form field Bmn through Hlmn = 3V[l-Ba/at]- 
Then the stability of the internal dimensions can be studied by the coupled 
three equations of motion involving gMN, (ft and H LMN instead of two (Eqs. 
13] and El). We expect an effective potential with a confining form so that the 
internal dimensional will oscillate and remain small. 

Classical supergravity holds when all radii are larger than the ten dimen- 
sional Planck length. So the equations we used are valid when the radii are 
either constant or growing with time. In this case we can safely neglect the 
massless excitations on the branes which will redshift quickly. We can also 




(55) 




(56) 
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neglect the brane antibrane annihilations because branes will be frozen as 
the transverse dimensions expand. However, in string theory, gravitational 
interaction is described not by metric alone but by the coupled system of 
metric and dilaton. The assumption that dilaton is constant may be incon- 
sistent with the cosmological equations describing the evolution. Also the 
production of stringy objects in the general time dependent background is 
complicated since the back reaction of such strings or branes is likely to be 
important. Further studies on the late-time behavior of brane gas cosmology 
with interacting branes together with dilaton field are needed. 
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